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Abstract. We consider an inverse problem for a hyperbolic partial differential 
equation on a compact Riemannian manifold. Assuming that Ti and T2 are two 
disjoint open subsets of the boundary of the manifold we define the restricted Dirichlet- 
to-Neumann operator Ar 1 ,r 2 - This operator corresponds the boundary measurements 
when we have smooth sources supported on i\ and the fields produced by these sources 
are observed on T2- We show that when Ti and Y2 are disjoint but their closures 
intersect at least at one point, then the restricted Dirichlet-to-Neumann operator 
Ar 1 .r 2 determines the Riemannian manifold and the metric on it up to an isometry. In 
the Euclidian space, the result yields that an anisotropic wave speed inside a compact 
body is determined, up to a natural coordinate transformations, by measurements 
on the boundary of the body even when wave sources are kept away from receivers. 
Moreover, we show that if we have three arbitrary non-empty open subsets Ti,^, 
and r,3 of the boundary, then the restricted Dirichlet-to-Neumann operators Ar-.,r fc for 
1 < j < k < 3 determine the Riemannian manifold to an isometry. Similar result is 
proven also for the finite-time boundary measurements when the hyperbolic equation 
satisfies an exact controllability condition. 
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1. Introduction and main results 



Let M be a compact and connected C^-smooth manifold of dimension n and let g be a 
C°°-smooth Riemannian metric on M. Let q be a real-valued C^-smooth function on 
M, and denote by A g the Laplace-Beltrami operator on M. We consider a hyperbolic 
inverse problem corresponding to the 2nd order elliptic operator 

a(x, D) := —A g + q(x). 
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In local coordinates g is a positive-definite C°°-smooth matrix {gjk{x))^ k= i with 
j,k= 



the inverse (fi ,J,fc (x))™ fc=1 and 



a(x,D)u=-\g\-W £ h {^^wA + ™ W 

j,k=l ^ ' 
where \g\ := det(^-fc). Hence our results cover the setting, where M C R" is an open 
domain with smooth boundary and a(x, D) is an elliptic operator of the form (pQ). 

Let H S (M) be the Sobolev space of s G N times weakly differentiable functions on 
M, and let Hq(M) be the H 1 (M) closure of Cq°(M), the space of smooth compactly 
supported functions. The operator 

Au(x) := a(x, D)u, D(A) := H 2 (M) n H^(M) (2) 

is self-adjoint in L 2 (M) = L 2 (M,dV g ), where dV g is the Riemannian volume measure. 
In local coordinates dV g = {g^^dx. 

Denote by v'{x,t) = v(x,t) the solution of the initial boundary value problem 

d 2 t v + a{x,D)v = inMx(0,oo), (3) 

v \dMx(p,oo) — /) 

v\ t =o = d t v\ t =o = 0, 

for / G C^°(dM x (0, oo)), and define the hyperbolic Dirichlet-to-Neumann operator 

A : C™(dM x (0, oo)) C°°(dM x (0, oo)), Af := d u v f \ dMxR+ , 

where d u is the normal derivative on dM. In local coordinates the exterior conormal v 
is the covector (y\, . . . , u n ) with 

n 

u j (x)g jk (x)u k (x) = 1, x G dM, 

j,k=i 



and 



n 8 



Denote by r2 the restriction of the Dirichlet-to-Neumann operator 

ALr 2 : C °°(ri x (0,T)) ^ C°°(r 2 x (0,T)), 

where ri,T 2 C dM are open. Furthermore, denote Ar lt r 2 ; — A-r\,r 2 - 

It is well known, that the map A determine the manifold (M, g) up to an isometry 
[5]. This is also true for the restriction Af r when V is nonempty and T is sufficiently 
large [23] . 

In many applications observations of physical fields can not be done on the same 
locations where the sources of the fields are. For instance, in imaging in Earth Sciences, 
elastic or acoustic fields are often implemented using explosions [101 [38] ■ In such a case 
observation devices need to be far away from the sources. 

Similarly, in electromagnetic imaging, it is technically difficult to use electrodes 
at the same time as sources and for making observations. These are typical examples 
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of cases where the observation devices and the sources of the fields are supported on 
disjoint sets. 

In this paper we show, that for certain collections of pairs (Ti,]^) of open and 
disjoint subsets of dM, the operators Ap r2 determine the manifold (M, g) up to an 
isometry. 

Theorem 1. Let r 1 ,r 2 ,S C dM be open sets such that Y\,Y2 C S and Y\ fl Y2 7^ 0. 
Then S ; given as a smooth manifold, and the operator h.T 1 ,T 2 determine the manifold 
(M, g) up to an isometry. 

Theorem 2. Let Y\, r 2 ,r 3 C dM be open and nonempty. Then the smooth manifolds 
T p , p = 1, 2, 3, and the operators 

Ari,r 2 > A rii r 3 , A r2i r 3 

determine the manifold (M, g) up to an isometry. 

For measurements on a finite time interval, we prove a theorem similar to Theorem 
[2] under an additional controllability assumption: 

(A) For any w G L 2 (M) there is a boundary value / G L 2 (dM x (0, 00)) satisfying 
v f (T/2)=w, supp(/)cr 3 x(0,T), 
where is the solution of the equation (j3J) and T 3 C dM. 

Let us comment the controllability assumption (A) when M is embedded in IR n . 
Property (A) follows from the geometric control condition of Bardos, Lebeau and Rauch 
[2], which yields exact controllability of the wave equation. Property (A) follows also 
from existence of a strictly convex function h on M with respect to the Riemannian 
metric g. 

Suppose that h G C 2 (M) is strictly convex and that p > is a lower bound for the 
Hessian of h in the Riemannian metric g, that is 

D 2 h(X, X) > p\X\ g , X G T X M, x G M. 

By j3U, (A) holds if 
4 

T > — max \Vgh\ g , sup V g h(x) ■ p(x) < 0, 

P xeM xeV 3 

where V 9 and | • \ g are the gradient and length with respect to the Riemannian metric 
c?, 17 is the Euclidean unit outward normal to DM C M n and V g h ■ V is the Euclidean 
inner product. We refer to [31] for examples of Riemannian manifolds (M,g) having a 
strictly convex function h. 

Theorem 3. Let r 1 ,r2,r 3 C dM be open and nonempty. If the controllability 
assumption (A) holds for T 3 and T > 0, and 

4d(x,y)<T, x G T p , y G M, p = 1,2, 3, 

then the Riemannian manifolds (T p ,g\^ ), p = 1,2,3, and the operators 

A T A T A T 
iV ri,r 2 > iV ri,r 3 ' yi r 2 ,r 3 

determine the manifold (M, g) up to an isometry. 
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Figure 1. On left, the measurements Ari,r a , ^-r u r 3 , A-r 2 ,r 3 are shown as arrows 
pointing from the support of sources to the support of observations. Using Lemma 
[T] we can change the direction of any arrow in the picture on left. Hence also the 
measurements shown on right are covered by Theorems [2] and El 

The proofs of these theorems consist of showing that the data determine, up to a 
gauge transformation, the boundary spectral data of the operator A on a part of the 
boundary. The manifold is then determined up to an isometry, as can be seen using the 
boundary control method [3j [5j |2U [25] . 

Notice also, that the operator Af r2 determines the operator Af Fl by a time 
reversal argument. 

Lemma 1. Let Ti,]^ C dM be open and nonempty. Let T > 0, and define the time 
reversal operator Rf(x,t) := f(x,T-t). Iff E C^(T 1 x (0, T)) andh e C °°(r 2 x (0, T)), 
then 

(/'^■r 2 ,ri^)L 2 (9Mx(0,T)) = (-RAp l5 p 2 i?/, 1l)L2(dMx(0,T))- 

Hence (Tj,g\-p.), j = 1,2, given as Riemannian manifolds, and the operator A r £ iT2 
determine the operator Af 2 Tl = (i?Ap x r2 i2)*. 

This result is relatively well known, see e.g. [E1II2], but for the sake of completeness, 
we will give a proof in the appendix. 

Let us review previous results on the topic. The inverse problem for isotropic wave 
equation on a compact manifold with measurements on the whole boundary was solved 
by Belishev and Kurylev [5]. This was based on the boundary control method originally 
developed in [3] for wave equation on a bounded domain of M. n . The inverse problems 
for more general hyperbolic equations on a compact Riemannian manifold with sources 
and observations on the same open subset T of the boundary has been studied by 
Katchalov and Kurylev [2l] , see also [29] . Similar problem has recently been studied for 
non-compact manifolds in [23| [26] . 

Inverse problems for elliptic equations with data on a part of the boundary have 
been studied intensively as they are the natural generalization of the Calderon's inverse 
problem for the conductivity equation [ID] . 
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When measurements are given on the whole boundary, the inverse problem for 
Schrodinger equation on a bounded domain of K n , n > 3, and hence for isotropic 
conductivity equation, was solved by Sylvester and Uhlmann in [39] . The corresponding 
two dimensional problem for isotropic conductivity equation was solved first by Nachman 
in [36] for C 2 conductivities, and for L°° conductivities, for which Calderon's inverse 
problem was originally posed, by Astala and Paivarinta in [TJ. Recently, also the 
inverse problem for Schrodinger equation on a bounded domain of dimension two with 
measurements on the whole boundary was solved by Bukgheim in [8] . The corresponding 
problem on a compact Riemannian surface was later solved in |17j . 

The inverse problem for Schrodinger equation on a bounded domain of MJ 1 , n > 3, 
with observations on an open subset T of the boundary was solved in [28]. The inverse 
problem for Schrodinger equation on a bounded domain of M. n , n = 2, with sources 
and observations on the same open subset T of the boundary was solved by Imanuvilov, 
Uhlmann and Yamamoto in [21]. The corresponding problem on a compact Riemannian 
surface was later solved in [16J. For related results with measurements on a part of the 
boundary , see US [22] . 

The inverse problem for the Laplace-Beltrami operator A g on a compact 
Riemannian manifold with sources and observations on the same open subset T of the 
boundary has been studied on analytic Riemannian manifolds of dimension n > 3 in 
[3U [32], and on Riemannian surfaces in [33J, see also [TBI IT9] . The inverse problem 
for the Laplace-Beltrami operator in dimensions n > 3 is open in general, even when 
measurements are given on the whole boundary. For positive results under certain 
geometrical conditions see |13j . 

2. Spectral analysis of the data 

Denote by (Xj)jeN the increasing sequence of distinct eigenvalues of the operator A and 
let (0fc)fceN be an orthonormal basis of real- valued C^-smooth eigenfunctions. Moreover, 
let be a partition of N such that ((pkjkeij is a basis for the space of eigenfunctions 

corresponding Xj. 

Let / G C C °°(<9M x (0, oo)), and consider Af also as a function in C°°{dM x R) by 
defining Af(-,t) = for t < 0. There is a constant C > such that for x G dM, the 
Fourier transform J r t ^ T Af{x) is an analytic function of r when Imr < — C. It is known 
(see e.g. [25]), that J r ^ r A/(a;) extends to a meromorphic function of r G C, and that 
it may have poles only at points y/Xj. Moreover, the residues at these points are 






res J r t ^ T Af(x) = S~] (/(•, \/X~j),d u <p k ) d„(p k (x), 



where f(x,r) = (J-t-yrf)(x, r) and dS g is the Riemannian surface measure. 
If j G N, a k are constants for k G Ij, and the linear combination 
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vanish on a nonempty open subset of dM, then = for all k G Ij by unique 
continuation, see e.g. [35]. Hence for an open nonempty set T C dM and j G N, 
the functions {d u (pk\r)k&ij are linearly independent. 

Let r!,r 2 C dM be open and nonempty. By linear independence and smoothness 
of d u (f>k, there are / G C£°(Ti) and x G T 2 such that 

fee/,- 

Moreover, for fixed r G C, the map / H> /(-,r) from C^fTi x (0, oo)) to C£°(Ti) is 
surjective. 

Hence the operator Ar l5 r 2 determines the eigenvalues Xj and the operators 

W 2 ;i = C c °°(ri x (0, oo)) -> C°°(r 2 x (0, oo)), (4) 

Lri,T 2 ;jf := S ^{fidv { Pk)L 2 {dM4S g )du¥k\T2- 

fee/,- 

3. Inverse problem with disjoint sources and observations 

In this section we prove Theorem [TJ 

Lemma 2. Suppose that f,hE C°°(R) are suc/j i/iai 

^(/»(*)/0r)/(y))|^ >lM , = #W*)/(*)/(v))I«m>,,m> 



for all j, fc G N. Then d j f(0) = for all j G N := {0, 1,2,.. .} or <9 fc /i(0) = /or a// 
positive k G N. 

Proof. Assume that the claim is not valid. Then there exist j G N and G N\ {0} such 
that <9 J /(0) 7^ and d k h(0) ^ 0. Let us next consider the smallest such integers j and 
k. 

By Leibniz's formula 

= dl +k dl(h(x)f(x)f(y))\ x=0 , y=0 - didi +k (h(x)f(x)f(y))\ x=0 , y=0 
j+k 



z=o 

.7 



m=0 



(^)^ m ^(o)^ m /(o)^' + V(o) 



where 



Si+[ J k )d k h(0)&>f(0)&>f(0) + S 2 - S 3 , 



Si := J2y + i ) d l h(0)& +k - l f(0)&f(0), 

s * := E C 1^)^(0)^7(0)^7(0), 
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S 3 := [ jd m h(0)d j - m f(0)d J+k f(0), 

and the terms with indices I = and m = have cancelled each other out. 

As k is the smallest positive integer such that d k h(0) ^ 0, we have d l h(0) = in 
the sum Si, and so Si = 0. As j is the smallest integer such that d^f(0) ^ 0, we have 
d j+k ~ l f(0) = in the sum S 2 and d^ m f(0) = in the sum S3, thus S 2 = S 3 = 0. 

Hence d /i(0)(c^/(0)) 2 = 0, which is a contradiction with the assumption that 
d j f(0) ^ and d k h(0) ^ 0. This proves the claim. □ 

In the proof of the next lemma we use the equation 

= E 4^/(*"K. ( 5 ) 

z — ' a! 

|a|=j 

where / G C°°(R"), t G R, u G R n and j G N, 

Lemma 3. Suppose that f,h& C°°(M. n ) are snc/j i/jai 

^(^/^/(y))^^ = d^(h(x)f(x)f(y))\ x=0 , y=0 

for all multi-indices a, /? G N n . JTien d a f(0) = /or a// multi-indices a G N n or 
d^h(0) = /or a// nonzero multi-indices (3 G N n . 

Proo/. Let j, fc G N and i> G R n . By © 

##(M^)/fa)/(™))U,^> - a t fe ^(/ i (tr)/(tr)/(^))| t=0iS=0 

= E E ^| U<V {P{hf)(p)dPf(p) - dP(hf)(0)d a f(0)) = 0. 

Hence 

^^(/ i (tr)/(tr)/( S r))| t=0 ,, =0 = #^M*«)/fa)/(H)lt=o,.=o 

for all j,k G N and r G R". 
Define the sets 

F := {r G R" : dP t f{tv)\ t=Q = for all j G N}, 

# := {v G R" : d k h(tv)\ t=0 = for all positive k G N}. 

The sets F and are closed by smoothness of / and h, respectively. Lemma [2] gives 
that F U H = R™. If F ^ R n , then R™ \ F is open, nonempty and contained in H. Thus 
F or H contains an open nonempty subset. 

Suppose that U C F is open and nonempty. Let j G N, and define the polynomial 

P(«) == E 4ff7(°)«°- 
z — ' a! 

|a|=J 

By (JSJ), p(r) = &[ f(tv)\ t =o, and p vanish in U. Using unique continuation for real 
analytic functions we see that p = in R™, and so the coefficients of p vanish. As j can 
be chosen freely, d a f(0) = for all multi-indices a. 
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Similarly, if there exists an open and nonempty V C H, then d a h(0) = for all 
nonzero multi- indices a. □ 

Remark 1. Let U be a -smooth manifold of dimension n, f G C°°(U) and p G U . 
If d a f(0) = for all multi-indices a G N n in some local coordinates taking p to 0, then 
d a f(0) = for all multi-indices a G N n in all local coordinates taking p to 0. 

Lemma 4. Let <fi be an eigenfunction of the operator A corresponding to an eigenvalue 
X, and let po G dM. Then in any local coordinates of dM taking po to 0, there is a 
multi-index a G N n_1 such that d a d u (j)(0) ^ 0. 

Proof. Assume that the claim is not valid. Then d a d l ,<f)(0) = for all a G N n_1 in some 
local coordinates of dM taking p to 0. 

Consider boundary normal coordinates of M taking po to 0. We may suppose 
that the coordinates map a small neigborhood V of p^ onto -8(0, e) x [0,e), where 
-8(0, e) C IR n_1 is a ball of radius e > centered at the origin. Then these coordinates 
take a boundary point p' G dM fl V to a point (x',0) G -8(0, e) x {0}, where 
x' = (x 1 , . . . , x n ~ l ). 

The special property of the boundary normal coordinates is, that a point p G MC\V 
has coordinates (a/, x n ) G -8(0, e) x [0, e), where x n = d(p, dM) and x' are the coordinates 
of the unique boundary point p' G dM such that d(p,p') = d(p, dM). 

Moreover, in the coordinates (x',x n ) the equation 

(-A g + q)4> = \<f>, 

has the form 

n— 1 n 

- d 2 ^ - 9 3k d x3 d xk ct> + £ a?d x ^ + aV = A0, (6) 
j,fc=i j=i 

for some a j G C°°(5(0, e) x [0, e)), j — 0, . . . , n, see e.g. [II] . 

Let us show, that = 0. Let b G N and a G N n_1 . By applying the operator d x ,d x „ 
on the both sides of flUD, we get 

n— 1 n 

«9^+ 2 0(O) = ^,^„(- ^ j^d^ + ^ + aV - A0)| x=o . 

j,fc=i j=i 

The right hand side of this equation is a linear combination of functions 

dgdht, a'eN,b'<b + l 
at the point x — 0. 

The equations c%</>(0) = hold for all a G N n_1 by the boundary condition 
4>\dM = 0. Furthermore, we have by Remark [fl that d x ,d x n(f)(0) = for all a G N n ~ x . 
Using induction we see that d x ,d x n(j)(0) = for all a G N n_1 and 6gN. 

Define the odd and even reflection operators 

R f(x',x n ) := (sign x")/(x', R e f{x',x n ) := fix', \x n \), 
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where / G C(B(O z e) x [0,e)), x' G B(0,e) and x n G (-e,e). 

Also, define := R a (f), g jk := R e g jk , a n := R Q a n and a? := R e a? for j = 0, . . . , n — 1. 
Denote ?7 := 5(0, e) x (-e,e). 

As 0|a;n=o = 0, we see that G H 2 (U) and 

d%d b xn 0(x', x n ) = (sign x n ) fc+1 cr,c£„0(x', |x n |), |a| + 6 < 2. 
Moreover, <p' fe is Lipschitz continuous in U, a J G L°°(U) for j = 0, . . . , n, and 

n— 1 n 

-dln(j) - 9 Jkd xid x k<f) + ^2a J d x j(f) + a°0 = A0, 

j,k=l j=l 

where the both sides are considered as functions in L 2 (U). Hence for some constant 
C > 

n-1 

j,k=l |q|+6<1 

Since G C°°(B(0,e) x [0,e)) vanishes up to arbitrary degree in origin, Taylor's 
formula gives for any m G N a constant C m > such that 

r r r 



|0(:r)|WW< / / C m r m dx n dx', asr->0. 

S(0,r) ./0 JB(0,r) JO 

Hence for any m G N, there is a constant C' m > such that 

\6{x)\ 2 dx n dx' < CLr m , asr^O. 



m i 

B(0,r) J-r 

By Hormander's strong unique continuation result [20J this yields, that = in U. 
In particular, = around some point q G M. As M is connected, unique continuation 
gives that = in M. This is a contradiction with the assumption that is an 
eigenfunction, and the claim is proved. □ 

Remark 2. Let X and Y be Hilbert spaces, and let Ux, ■ ■ ■ , Wjv G X and v i, . . . , vn G Y 
be linearly independent. Suppose that D C X is a dense subspace, and define 

N 

L:D^Y, Lf:=J2(f,u k ) x v k . 

k=l 

Then L determines the unique bounded extension L : X — » Y , and its adjoint 
L* : Y — ^ X . Hence, L determines the spaces 

span(t>i, . . . , vn) = L(X), span(«i, . . . , un) = L*(Y). 

Theorem 4. Let Tx,V 2 ,E C dM be open, T 1 ,T 2 C E, and fx n f 2 ^ 0. Then the 
smooth manifold £ and the collection 

{U,-/., ,r,: ; ) I JGN} (7) 
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determine boundary spectral data up to a constant gauge transformation on T 2 ■ That is, 
one can find a collection 

{(A is (^Vfclrjfce/,-) I 3 e N}, (8) 

where for an unknown constant C > not depending on j or k, (Cif)k)k€l- ^ s an 
orthonormal basis of eigenfunctions in L 2 (M) corresponding the eigenvalue Xj. 

Proof. Choose a smooth positive measure d/i on E. Then there is a positive function 
rj G C°°(S) such that rjd/i = dS g \^. As r\ > 0, the functions {T)d v <j>k\ri)keij are linearly 
independent for all j G N. 

For all j G N, denote Lj := L rit Y 2 -j, and define 

lj(x, y) := 2j r){x)d v <j)k{x)d v <j>k{y), x, y G dM, 

E) : = span(?7^0 fc |r 1 )fcG/ J , Ej ■= span (<9„0 fc |r 2 

Note, that for the smallest eigenvalue Ao, the space of eigenfunctions is one 
dimensional (see e.g. p3| Thm 6.5.2]), and so 

l {x,y) = ri(x)d u (f> (x)d v (f)o{y)- 

Consider a positive function rj G C°°(S) and real- valued functions G C°°(E), 
G N, such that the following three conditions hold: 

(Al) If xq G Ti n T2, then in local coordinates of £ taking xq to 

™ r& ( k(x,y) \ , _ ^™ / /o(a:,y) \ , 

d * d * \~WTJ lx=0 ' y=0 - dxdy VWV) lx=0 ' y=0 

for all multi-indices a, f3 G N n_1 . 
(A2) span^ejfclrjjfceii = #j an d span(e fc |r 2 )fce/ j = E) for all j G N. 
(A3) L, = Lj for all j G N, where 

L jf(y) '■= "£2(f>vek)L2(dM,d») e k(y), f e Co°(ri), y g r 2 . 

fee/. 

Such functions 77 and e k exist. For example, rj = rj and e& = d u (pk\T, satisfy the 
conditions. 

Next we show the following two statements. 

(i) We can verify using the data ([7]), whether any given functions rj G C°°(S) and 
e k G C°°(E), keN, satisfy the conditions (Al), (A2) and (A3). 

(ii) There is an orthonormal basis of eigenfunctions (ipk)kem of operator A and a 
constant C > 0, not depending on k, such that 

e fe |r 2 = Cd u ip k \r 2 , k G N. 
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If (i) holds, then the data fl7|) determine the nonempty collection 

{(efclrJfceN | (A1)-(A3) hold with a positive rj E C°°(£)}, 

and if (ii) holds, then any element from this collection determines a collection of type 
OH]). So the claim of the theorem is proved after proving (i) and (ii). 

Let us show the claim (i). Clearly, the condition (A3) can be verified using the data 
©. As 

L 2 (dM,dfj,), f G C™(Ti),y E r 2 , 

by varying / E C^°(r 1 ), we see that the map Lj determines the function lj\r 1 xr 2 - 
Let xq E Ti PI T 2 . Given Zo|rixr 2 an d rj, it is possible to compute 



in any local coordinates £ taking x to 0. By smoothness of t], rj and <9j,0cb these 
functions are known also at (x, y) = (0,0). Hence the condition (Al) can be verified 
using the data (j7j). 

Taking X = L 2 (T u d/j>) } Y = L 2 (T 2} dfi) and D =JJ™(T 1 ) in the formulation of 
Remark [21 we see that the map Lj determines the spaces Ej and Ej. Hence the condition 
(A2) can be verified using the data (JTj), and the claim (i) is proved. 

Let us show the claim (ii). Let xq E rifl^. LemmaHlgives that, in local coordinates 
of DM taking xq to 0, there is a multi-index a E N n_1 such that d a d u (po(0) ^ 0. Hence 
the condition (Al) and Lemma |3] imply, that <9^(?7r/ _1 )(0) = for all nonzero multi- 
indices E N n -\ 

Fix j E N and, to simplify the notation, drop the subindices j from now on. By 
the condition (A2) 

r l e i\r 1 = / ] aiky9 v (j)k\rn e z |r 2 = 5j&a^Jfc|r 2 I 6 /, 
kei kei 

for some constant matrices A := (aik)i y kel an d B := (bik)i,kel- 
Fix xq E Ti n r 2 , let / E I and define the function 



We have seen that, in local coordinates of dM taking xq to 0, the equation (^(7777 x )(0) 
holds for all nonzero multi-indices (3 E N n_1 . Hence for any multi- index a E N n ~ x 

= d a e t (0) -d a ei(0) 



aiJL^-rd^kix) 



\kei 

£ - ^) 9°^fc(0) = <m0(O). 



x=0 V fc G 7 



y=0 
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By Lemma HI the coefficients 

aik^T^r-kk, k,lel 

7/(0) 

vanish, and so 77(0)77(0) _1 A = B. 
Moreover 

(h, Lf) L 2(QM,dS g ) 

= /jifi V e k)L 2 (dM,d^)(h, e k ) L 2( dMjdSg ) 
kei 



m) L 2 (dM,dS g ) 

kei lei mei 

m) L 2 (dM,dS g )i 

l,mel \kel J 

for all / G C™(Ti) and h G C£°(r 2 ). On the other hand, the condition (A3) gives 
{h, Lf) L 2(QM,dS g ) = (h, Lf) L 2( dM>dSg ) 

L 2 (dM,dS g ){h, d 1/ (j)k)L 2 (dM,dSg)^ 

kei 

for all / G Cg°(ri) and h G Cg°(r 2 ). 

Denote (-, ■) := (-, -)L 2 (dM,ds 3 )- By density of Cg°(r p ) in L 2 (T P , dS g ), p = 1, 2, 

5^ I 5^ aM&fcm ) U,dv<l>l)(KdAm) = ^2{f,d u <Pk){h,d u (f) k ), 
i,mei \kei J kei 

for all / G L 2 (r 1; dS g ) and /i G L 2 (r 2 ,dS 9 ). 

Let (fi)i e i be biorthogonal with (d„<f>i\ Vl ) leI in L 2 (ri, dS 9 ), and let(/i m ) TO&f be 
biorthogonal with (d v (j) m \r 2 )mei in L 2 (T 2 , dS g ), that is, 

(fi>, d u <pi) = 8m, (h m >, d u <fi m ) = 6 m > m , I', I, m, w! G I. 

Then 

O-kl'bkm' = I '^^Cbklbkm J {fl f , d v (j)j) (h m > , d v (j} m ) 

kei i,mei \kei J 

= y^Xfv, d„(j)k){h m >, d u <f>k) = <W, l', m ' e J. 
kei 

Denote c := r/(0)~ 1 77(0) > 0. We have shown, that / = A T B = cB T B. Hence the 
matrix y/cB is orthogonal. To conclude, we observe that 



r 2 = —j=dv y~] Vcb lk (j) k \ r2 , I G /, 
\/c z — ' 
v kei 



where (J2kei V^ik^kjiel is an orthonormal basis of eigenfunctions corresponding the 
eigenvalue A,-. □ 
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As discussed in the previous section, the operator Ar 1) r 2 determines the collection 
(J7|). So by previous theorem, if Ti fl T 2 7^ 0, then the operator Ap lj r 2 determines the 
collection (jSJ). The collection (JH]) determines the manifold up to an isometry by [25, 
Chapter 4.4]. This proves Theorem [TJ 

4. Inverse problem with observations far away from sources 

In this section we prove Theorems [2] and 

Proof of Theorem^ Denote L p ^ q = Lr pt r q -,j- It is enough to show, that the collection 

{(A„L^ 2 ,L]- 2 ,L^ 3 ) I jeN} (9) 

determines a collection of type (jHJ). 

Denote £ = Fx U T 2 U T 3 , and choose a smooth positive measure dfi on E. There is 
a positive function 77 G C°°(E) such that 77c?// = gLSJe. Define for all j G N 

£j := span^d^fclrjfcg/,, p = 1, 2, 

E] : = span(^0 fc |r q )fcG/ J) g = 2,3. 

Choose a positive function 77 G C°°(S) and real- valued functions G C°°(S), 
fc G N, such that the following two conditions hold: 

(Bl) for all j G N 

span^efclrjfce/, = Ej, V = 1, 2, 
span(e fe |r g )fe e / 3 = J5?, 9 = 2,3. 

(B2) Z^ 9 = Ly q for all j G N and (p, q) = (1, 2), (1, 3), (2, 3), where 

Again, such functions 77 and exist, as rj = 77 and e k = c^^ls satisfy the conditions 
(Bl) and (B2). It is enough to show the following two statements: 

(i) We can verify using the data (02), whether any given functions rj G C°°(E) and 
e k G C°°(E), fc G N, satisfy the conditions (Bl) and (B2). 

(ii) There is an orthonormal basis of eigenfunctions (ipk)kem of operator A and a 
constant C > 0, not depending on fc, such that 

e k \r 2 = Cd u ^ k \ r2 , keN. 

Analogously with the proof of Theorem HI the maps L^ 2 , L]~^ 3 and E 2 ^ 3 determine 
the spaces Ej, E?, E 3 and E 2 . Hence the claim (i) is proved. 

Let us show the claim (ii). Fix j 6 N and, to simplify the notation, drop the 
subindices j from now on. The condition (Bl) gives, that for all / G / 

kei fee/ 

e '|r 2 = ^J^fe^0fe|r 2 ) e dr 3 = ^^Qfect0fe|r 3 , 
fee/ fee/ 
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for some constant matrices 

A '■= { a ik)i,k&ii B := (bik)i,kei, B := (bi k )i tk€l , C := (ci^i^ei- 
For the smallest eigenvalue Ao, the space of eigenfunctions is one dimensional, and 

so 

(ri(y)~ \ 
x-r-^oo - b 00 J d u <po{y) 

for all y G T 2 . By Lemma HI the set 

N := {y G T 2 : d u Mv) = 0} 



does not contain a nonempty open set of dM. Hence T 2 \ = T 2 . Moreover, 
7777" 1 £>oo — ^00 = in T 2 \ N and by continuity also in the whole set T 2 . 

Denote by c the constant r]~ 1 rj\r 2 > 0- As (d u (j) k \r 2 )kei are linearly independent, 
B = cB for all j G N. Moreover, we may use the condition (B2) as we used the 
corresponding condition in the proof of Theorem HJ and get 

A T B = I, A T C = I, B T C = I. 

Hence B = C and cB T B = I. 

To conclude, we observe that 



e i\r 2 = —j=dv Vcbik<t>k r 2 > 1 e I, 
v kei 



where {^2 keI \/cbik<pk)iei is an orthonormal basis of eigenfunctions corresponding the 
eigenvalue Xj. □ 

We prove Theorem[3]by reduction to Theorem [2] using a time continuation argument 
similar to 



Lemma 5. Suppose that Ti,r 2 C dM are open and nonempty. Denote 

T* :=2max{d{x,y) : x G Ti, y G M}. 

If T* < t < T, then the smooth manifolds Ti, T 2 , the operator Af r2 and the inner 
products 

(u f {t ),u h (t )) L 2 iM) , f,he C™^ x (0,T)) (10) 
determine the operator Ap+p 2 for 5 < to — T* . 

Proof. Denote by Y s the time delay operator 

Y s f(-,t):=f(;t-s), t,seR. 

As the coefficients of the wave equation (j3J) are time-independent, v Ys *(x,t) = 
(Y s vf)(x,t) and (AY s f)(x,t) = (Y s Af)(x,t). 

Let / G C °°(ri x (0,T + 5) and choose h G C^(T 1 x (0,t )) and h! G C£°( r i x 
(5,T + 5)) such that f = h + h'. 
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Let e G (0,5). As swpp(Y_gh') C r x x (0, T), the operator ^f 1 r 2 determine the 
function 

Ah'(; T + e) = (y_iAA')(-, T - (5 - e)) = (A^Y.gh')^ T — (5 — e)), 

in T 2 . Therefore, it is enough to show that the given data determine also 

Ah(x,T + e), x G r 2 . 

Consider a sequence {hj)'JL l C C^°(ri x (0, to — ^)) satisfying the following two 
conditions. 

(CI) lim^^^^o) = « h (*o) in H\M), 
(C2) lim^oo^t;^^) = d t v h (t ) in L 2 (M). 

Such a sequence exists, since to > T*, and thus we see exactly as in (25J Thm. 4.28] 
that the set 

{(^(t - 5),<V(t ~S)):fe CfFx x (0,t - 5))} 

is dense in H%(M) x L 2 (M). 

Let us prove, that (CI) is equivalent with 

(CI') For all c>0 

lim ((-AgWj + qw j ,Wj) L 2r M) + c(w h Wj) L 2 {M) ) = 0, 
where Wj := v Yshl (to) — v h (t ). 

As supp(F^/ij) C Ti x (0, to) and supp(/i) C Ti x (0, to), we have that Wj\dM = 0. 
Hence 

-(AgWj,Wj) L 2( M ) = (dwj,dwj) L 2( M }, 
where d is the exterior derivative on M. If (CI) holds, then 

K-AgWj + qWj,Wj)unM) + c(w j} Wj)i^^\ 
< \\dWj\\ L 2( M -) + ||? + c|| LO o( M ) ll^ill^^) 
—> 0, as j — > oo. 

For large enough c > there is a constant c > such that q + c > c . Hence if (CI') 
holds, then 

IM W jIIl2(M) + C IK'IIl2(M) - (~ A g W j + QWj,Wj) L 2 {M) + c(Wj,Wj) L 2^ M) -»■ 0, 

as j — > oo, and (CI) holds. Therefore (CI) and (CI') are equivalent. 
Next we observe that 

_ a 2 (u y s (r^,-/ l ) (io)) ^,- /l(to))i2(M) | s=o 

= (-<9 f 2 u; j ,U7 i ) L 2 (M) = ((-A g + q)w j ,w j ) L 2 (M) . 

Hence the condition (CI') can be verified for given functions (hj)^L 1 and h using the 
inner products (llOj) . 
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Similarly, 

d Sl d S2 (v Y ^ Y ^- h \t ),v Y ^ h ^ 

and condition (C2) can be verified for given functions {hj)JL 1 and h using the inner 
products (jTOl . 

As v Yshi — v h = on dM x [t ,oo), conditions (CI) and (C2) together with the 
continuous dependency of the solution of the wave equation on the initial data, see e.g. 
[251 Thm. 2.30], give 

lim d v (v Ysh i - v h ) = 0, in L 2 (dM x (t , oo)). 

We have seen that, the inner products (TIP"]) determine for any h G C£°(Ti x (0, to)) 
the nonempty set 

{{hj)f =l C C?(T X x (0,t - 5)) I (CI) and (C2) hold}, 

and that any sequence in this set satisfies 

Ah(x,T + t)= ]im(AY s hj)(x,T + t)= lim r hj(x,T — (8 — t)) 

in L 2 (r 2 x (T,T + 5)). As A/i G C°°(r 2 x (0,T + 5)), the inner products (JTHJ and the 
operator Ap p 2 determine Ah(x, T + e) pointwise for x G T 2 and e G (0, 5). □ 

Next we prove the last of the three main theorems formulated in the introduction. 

Proof of Theorem^ By Lemma [1] the operators Ap p 2 ,Ap r , Af r determine the 
operators 

A^p q , p, g = l,2, 3, p^g, (11) 
We use the time delay operator Y s defined in the proof of Lemma [51 Define 

B[f, h] := I [ (d u v^ - v f d^A dS g dt, f,he C™{dM x (0, oo)), 

JO JdM ^ ' 

and let to := T/2. We recall Blagovestchenskii identity [251 Lem. 4.16], originating from 

®, 

1 f to 

(v f (to) 1 v h (t )) L 2 {M) = - / (sign s)B[Y to+s f,Y to _ s h]ds, 

1 J -t 

where /, h G C^°(dM x (0,T)). By this identity, the operators fllip determine the inner 
products 

(v f (t ),v h {t )) L 2 {M) , f G C^(T P x (0,T)), h G ^(r, x (0,T)), (12) 

for p,g = 1,2,3, p ^ g. 

Next we will show that the operators ( ITT]) determine the inner products (|T2|) also 
for p = g, p = 1, 2, 3. 

Let / G C °°(r 2 x (0, T)) and consider a sequence C C °°(r 3 x(0, T)) satisfying 

the following two conditions. 
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(Dl) For all h G C^(T 1 x (0,T)) 

lim (v f (t ) - v^(t ),v h (t )) L 2 [M) = 0. 

(D2) The sequence {fj)°° =1 is bounded in L 2 {T 3 x (0,T)). 

By assumption (A), there is / G L 2 (T 3 x (0,T)) such that v f (t ) = v^(t ). Thus, 
there is a sequence (/ 3 -)£Li C Co°(r 3 x (0, T)) such that 

lim/, = /, inL 2 (r 3 x (0,T)). 

By |2SJ Lem. 2.42], v^(to) — > vf(to) in L 2 (M) as j — >■ oo. Hence a sequence satisfying 
the conditions (Dl) and (D2) exists. 
Let us next show that 



K(*o) 



lim (f /j (to),f / (^o))L2( M )- 



il 2 (m) /: 

As t > 2d(x,y) for all x G I\ and y G M, [251 Thm. 3.10] gives that the set 

{v^lheC^xfaT))} 
is dense in L 2 (M). Let e > and choose h G Co°(Ti x (0,T)) such that 



(13) 



\v f (t )-v h (t )\ 



L 2 (M) 



< e. 



By [25| Lem. 2.42] there is C > 0, and by the condition (Dl) there is J G N such 
that for j > J 

/ L 2 (M) 



< 



+ 

< C 

By the condition (D2) 



v f (t )-v^(t Q ),v f (to)-v h {t ) 
v f(t )-ifi{t ),v h {t ) 



L 2 (M) 



L 2 (M) 



f-fj 



L 2 (dMx(0,T)) 



e + e. 



sup 



f-fi 



< OO. 



L 2 (dMx(0,T)) 

Hence the equation (TlBl is valid. 

By |25l Thm. 3.10] the functions v f (t ), f G C °°(r 2 x (0,T)), are dense in L 2 (M). 
Hence 



A*o)|| L 2 (M) = SUp(^(t ),^ / (^o))L2(M) 



(14) 



where h G C^(T p x (0,T)), p = 1,3, and the supremum is taken over all / G 
Cnr 2 x (0,T)) such that |K(*o)|| i2(M) = 1. 
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The condition (Dl) can be verified for any / and {fj)JLi using the inner products 
(j!2p for p = 2,3, q = 1. Therefore, these inner products determine for any / G 
C^(T 2 x (0,T)) the nonempty set 

C C °°(r 3 x (0,T)) I (Dl), (D2) hold}. 

By equation ([TBI any sequence in this set together with inner products (!T2|) for p = 3 
and g = 2 determine ||i/(£o) || L 2( M )- 

As / G C^°(r 2 x (0, T)) can be chosen freely, the inner products ( !P2]) for p = 2,3, 
q = 1 and for p = 3, q = 2 together with polarization identity determine the inner 
products ( IT2|) for p = g = 2. 

The equation (fT4l) . polarization identity and the inner products ( {121 for p — 1,2, 3, 
g = 2 determine the inner products (JT2J) for p = g = 1, 3. 

Therefore, the operators ffTTj) determine the inner products 

(v f (t ),v h (t )) L 2 {M) , f,heC™(T p x(0,T)), p= 1,2,3. (15) 

Choose 5 G (0,to — ^*); where T* is defined as in Lemma [5l By Lemma [5] the 
operators (TTTj) and the inner products (fT5l) determine the operators 

A^ q , p,g = l,2,3, p^g. (16) 

Repeating this construction, we see that the operators 

A£3?, P, 9 = 1,2, 3, p^g, 

are determined for all m G N. The claim follows from Theorem [2j □ 

Appendix 

Next we prove the Lemma [TJ stating that the operator Ap Fa determines the operator 
7V r 2) ri- 

Proof of LemmaUl Define u := Rv R * , where v R * is the solution of the equation ([3]) with 
the boundary data Rf G Co°(ri x (0,T)). Then u(x,t) = v R f(x,T — t) satisfies the 
equation 

d 2 t u + a(x, D)u = 0, in M x (0, T), 

U\dMx(0,T) = f, 

u\ t =T = d t u\ t=T = 0. 
Integration by parts gives 

(/> Ap 2 ri /i) L 2( 9Mx ( r )) - (RAr i r2 Rf, /i)L 2 (9Mx(o,T)) 

= f [ (u{x,t)d v v h {x,t) -(d v v Rf )(x,T -t)v h (x,t))dS g (x)dt 

J0 JdM 

= / (u(x,t)A g v h (x,t) - (A g u)(x,t)v h (x,t)) dV g (x)dt 
Jo Jm 



Inverse problem for wave equation with sources and observations on disjoint sets 19 

r-T 



M 



(u{x,t)d*v h (x,t) - (d 2 t u)(x,t)v h (x,t)) dV g (x)dt 
(u(x,t)d t v h (x,t) - (d t u)(x,t)v h (x,t)) dVg(x) 



t=0 



since u\ t =T = d t u\ t= T = and v^\ t =o — d t v^\ t =o = 0. □ 

Acknowledgements. The authors were partly supported by Finnish Centre of 
Excellence in Inverse Problems Research, Academy of Finland COE 213476. LO was 
partly supported also by Finnish Graduate School in Computational Sciences. 

References 

[1] K. Astala, L. Paivarinta, Calderon's inverse conductivity problem in the plane. Ann. of Math. (2) 

163 (2006), no. 1, 265-299. 
[2] C. Bardos, G. Lebeau, J. Rauch, Sharp sufficient conditions for the observation, control and 

stability of waves from the boundary, SIAM J. Control Optim. 30, 1992, 1024-1065. 
[3] M. Belishev, An approach to multidimensional inverse problems for the wave equation. (Russian) 

Dokl. Akad. Nauk SSSR 297 (1987), no. 3, 524-527; translation in Soviet Math. Dokl. 36 (1988), 

no. 3, 481-484. 

[4] M. Belishev, The Calderon problem for two-dimensional manifolds by the BC-method. SIAM J. 

Math. Anal. 35 (2003), no. 1, 172-182. 
[5] M. Belishev, Y. Kurylev, To the reconstruction of a Riemannian manifold via its spectral data 

(BC-method), Comm. Partial Differential Equations 17, 1992, no. 5-6, 767-804. 
[6] K. Bingham, Y. Kurylev, M. Lassas, S. Siltanen, Iterative time reversal control for inverse 

problems. Inverse Problems and Imaging 2 (2008), 63-81. 
[7] A.S. Blagovestchenskii, The inverse problem of the theory of seismic wave propagation. (Russian) 

Probl. of Math. Phys., No. 1, pp. 6881. Izdat. Leningrad. Univ., Leningrad, 1966. 
[8] A. L. Bukhgeim, Recovering a potential from Cauchy data in the two-dimensional case. J. Inverse 

Ill-Posed Probl. 16 (2008), no. 1, 19-33. 
[9] A. L. Bukhgeim, G. Uhlmann, Recovering a potential from partial Cauchy data. Comm. Partial 

Differential Equations 27 (2002), no. 3-4, 653-668. 
[10] A. P. Calderon, On an inverse boundary value problem. Seminar on Numerical Analysis and its 

Applications to Continuum Physics (Rio de Janeiro, 1980), pp. 65-73, Soc. Brasil. Mat., Rio de 

Janeiro, 1980. 

[11] I. Chavcl, Riemannian geometry: A modern introduction, Cambridge Studies in Advanced 

Mathematics 98, Cambridge University Press, 2006, xvi+471 pp. 
[12] M. Dahl, A. Kirpichnikova, M. Lassas, Focusing waves in unknown media by modified time reversal 

iteration. SIAM J. Control Optim. 48 (2009), no. 2, 839-858. 
[13] D. Dos Santos Ferreira, C. Kenig, M. Salo, G. Uhlmann, Limiting Carleman weights and anisotropic 

inverse problems. Invent. Math. 178 (2009), no. 1, 119-171. 
[14] L. Evans, Partial differential equations. Graduate Studies in Mathematics 19, American 

Mathematical Society, 1998, xviii+662 pp. 
[15] A. Greenlcaf, G Uhlmann Local uniqueness for the Dirichlet-to-Neumann map via the two-plane 

transform Duke Math. J. 108 (2001), 599-617. 
[16] C. Guillarmou, L. Tzou, Calderon inverse Problem with partial data on Riemann Surfaces. 

larXrv:0908.1417V 2. 

[17] C. Guillarmou, L. Tzou, Calderon inverse Problem for the Schrodinger Operator on Riemann 
Surfaces. larXiv:0904.3804. 



Inverse problem for wave equation with sources and observations on disjoint sets 20 



G. Henkin, V. Michel, On the explicit reconstruction of a Riemann surface from its Dirichlet- 
Neumann operator. Geom. Funct. Anal. 17 (2007), no. 1, 116-155 

G. Henkin, V. Michel, Inverse conductivity problem on Riemann surfaces. J. Geom. Anal. 18 
(2008), no. 4, 1033-1052. 

L. Hormander, Uniqueness Theorems for Second Order Elliptic Differential Equations, Comm. in 

Partial Differential Equations, 8(1), 21-64, 1983. 
O. Imanuvilov, G. Uhlmann, M. Yamamoto, Global uniqueness from partial Cauchy data in two 

dimensions. larXiv:0810."2286V l. 
V. Isakov, On the Uniqueness in the Inverse Conductivity Problem with Local Data. Inverse 

Problems and Imaging, Vol. 1 (2007), no. 1, 95-105. 

H. Isozaki, Y. Kurylev, M. Lassas, Forward and Inverse scattering on manifolds with asymptotically 
cylindrical ends. To appear in Journal of Functional Analysis. doi:10.1016/j.jfa.2009.11.009 

A. Katchalov, Ya. Kurylev, Multidimensional inverse problem with incomplete boundary spectral 

data, Comm. Partial Differential Equations 23, 1998, no. 1-2, 55-95. 
A. Katchalov, Y. Kurylev, M. Lassas, Inverse Boundary Spectral Problems, Monographs and 

Surveys in Pure and Applied Mathematics 123, Chapman Hall/CRC-press, 2001, xi+290 pp. 
A. Katchalov, Y. Kurylev, M. Lassas, Energy measurements and equivalence of boundary data 

for inverse problems on non-compact manifolds. IMA volumes in Mathematics and Applications 

(Springer Verlag) "Geometric Methods in Inverse Problems and PDE Control" Ed. C. Croke, I. 

Lasiecka, G. Uhlmann, M. Vogelius, 2004, pp. 183-214. 
A. Katchalov, Y. Kurylev, M. Lassas, N. Mandache, Equivalence of time-domain inverse problems 

and boundary spectral problems. Inverse Problems 20 (2004), no. 2, 419-436. 
C. Kenig, J. Sjostrand, G. Uhlmann, The Calderon problem with partial data. Ann. of Math. (2) 

165 (2007), no. 2, 567-591. 
Kurylev, M. Lassas, Hyperbolic inverse problem with data on a part of the boundary. Differential 

equations and mathematical physics (1999), 259-272, AMS/IP Stud. Adv. Math., 16, Amer. 

Math. Soc, 2000. 

Kurylev, M. Lassas, Hyperbolic inverse boundary-value problem and time-continuation of the 
non-stationary Dirichlet-to- Neumann map. Proc. Roy. Soc. Edinburgh Sect. A 132 (2002), no. 
4, 931-949. 

I. Lasiecka, R. Triggiani, P. Yao, Exact controllability for second-order hyperbolic equations 
with variable coefficient-principal part and first-order terms. Proceedings of the Second World 
Congress of Nonlinear Analysts, Part 1 (Athens, 1996). Nonlinear Anal. 30 (1997), no. 1, 111— 
122. 

M. Lassas, M. Taylor, G. Uhlmann, The Dirichlet-to-Neumann map for complete Riemannian 

manifolds with boundary. Comm. Anal. Geom. 11 (2003), no. 2, 207-221. 
M. Lassas, G. Uhlmann, On determining a Riemannian manifold from the Dirichlet-to-Neumann 

map. Ann. Sci. cole Norm. Sup. (4) 34 (2001), no. 5, 771-787. 
J. Lee, G. Uhlmann, Determining anisotropic real-analytic conductivities by boundary 

measurements. Comm. Pure Appl. Math. 42 (1989), no. 8, 1097-1112. 
R. Leis, Initial-boundary value problems in mathematical physics. B. G. Teubner, Stuttgart; John 

Wiley k Sons, Ltd., Chichester, 1986. viii+266 pp. 
A. Nachman, Reconstructions from boundary measurements. Ann. of Math. (2) 128 (1988), no. 3, 

531-576. 

A. Nachman, J. Sylvester, G. Uhlmann, An n-dimensional Borg-Levinson theorem. Comm. Math. 

Phys. 115 (1988), no. 4, 595-605. 
Rakesh, A linearised inverse problem for the wave equation. Comm. Partial Differential Equations 

13 (1988), no. 5, 573-601. 
J. Sylvester, G. Uhlmann, A global uniqueness theorem for an inverse boundary value problem. 

Ann. of Math. (2) 125 (1987), no. 1, 153-169. 
W. Symes, The seismic reflection inverse problem. Inverse Problems 25 (2009), no. 12. 



